AIAA-83-1 691 

Interactions  of  Airfoils  with 
Gusts  and  Concentrated  Vortices 
in  Unsteady  Transonic  Flow 

5  W.  J.  McCroskey  and  P.M.  Goorjian 
C;  NASA  Ames  Research  Center 
Moffett  Field,  CA 


'  r  .^’.-4  —  '—'T  u- 

OCT  2  3  1383  .>  | 


AIAA 1 6th  Fluid  and 
Plasma  Dynamics  Conference 


& 


July  12-14, 1983 
Danvers,  Massachusetts 


Q3  1 0  1  ** 


184 


sa 


ft 


m 


INTERACTIONS  OF  AIRFOILS  WITH  GUSTS  AND  CONCENTRATED  VORTICES  IN  UNSTEADY  TRANSONIC  FLOW 

W.  J.  McCroskey*  and  P.  M.  Goorjian" 

NASA  Ames  Research  Center,  Moffett  Field,  California 


Abstract 

Unsteady  interactions  of  concentrated  vortices 
and  distributed  free-stream  gusts  with  a  stationary 
airfoil  have  been  analyzed  in  two-dimensional  tran¬ 
sonic  flow.  A  simple  method  of  introducing  such 
disturbances  has  been  implemented  numerically  in 
the  well-known  transonic  small-disturbance  code 
LTRAN2,  and  calculations  have  been  performed  for 
two  important  classes  of  current  aerodynamic  prob¬ 
lems.  The  first,  which  demonstrates  many  of  the 
essential  features  of  the  interactions  between 
helicopter  rotor  blades  and  their  trailing-vortex 
wakes,  is  that  of  a  discrete  potential  vortex  con¬ 
verting  past  an  airfoil.  The  second  is  the 
response  of  a  transonic  airfoil  to  a  transverse 
periodic  gust,  with  and  without  the  alleviation 
that  can  be  achieved  by  the  proper  active  control 
motion  of  a  trailing-edge  flap.  In  both  cases, 
unsteady  effects  are  found  to  play  important  roles 
in  the  shock-wave  motion,  in  the  overall  flow-field 
development,  and  consequently,  in  the  air  loads  on 
the  airfoil, 
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Nomenclature 


=  time  index 


q  =  nondimensional  velocity 

<JV  =  velocity  induced  by  the  vortex 

^  =  velocity  at  the  location  of  the  vortex, 

excluding  qv 

r0  =  initial  position  of  the  vortex 
t1  =  instantaneous  position  of  the  vortex 

t, t  =  dimensionless  and  dimensional  time,  respec¬ 

tively,  sec 

U„  =  free-stream  velocity,  m/sec 

u, v  =  dimensionless  velocities  in  the  x-  and 

y-dlrections,  respectively 

uv,vv  «  velocities  induced  by  the  vortex 

Vq  =  gust  velocity 

x,y  =  directional  coordinates 


amplitude  ratio,  defined  following  Eq.  (14)  x0,y0  =  initial  position  of  the  vortex 


vortex  core  radius 


(1  -  Ha,)€~l/3 

-(l  +  y)m£ 


lift  coefficient 


moment  coefficient 


pressure  coefficient 
airfoil  chord,  m 

difference  operator,  defined  following 
Eq.  (3) 

defined  following  Eq.  (3) 

indices  in  x-  and  y-directions,  respec¬ 
tively;  also  unit  vectors  in  x-  and 
y-directions 

reduced  frequency,  wc/u*, 

Mach  number 

exponent  in  Eq.  (1);  also  dummy  time  index 
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xl>Yl  =  Instantaneous  position  of  the  vortex 

a  =  angle  of  attack,  deg 

y  ■  rafi  of  specific  heats 

T  "  strength  of  the  vortex 

Ax  «  difference  operator,  defined  following 
Eq.  (3) 

6  -  difference  operators,  defined  following 

Eq.  (3) 

6p  ■  flap  deflection,  deg 

c  •  airfoil  thickness  ratio 

8,0  =  phase  angles,  defined  following  Eq.  (14) 

A  ■  sweep  angle,  Fig.  1 

*  disturbance  potentials,  with  and  without 
the  vortex,  respectively 

w  ■  frequency  of  oscillation,  rad/sec 


I.  Introduction 

The  unsteady  interaction  of  concentrated  vor¬ 
tices,  gusts,  and  related  vortical  disturbances 
with  lifting  surfaces  can  be  important  for  a 
variety  of  fluid-dynamic  devices,  such  as  maneu¬ 
vering  aircraft  and  missiles,  helicopter  rotor 
blades,  and  turbomachinery .  Frequently,  these 
interactions  alter  the  flow  field  significantly, 
with  consequent  changes  in  aerodynamic  loading, 
aeroelastic  behavior,  and  acoustic  characteristics. 
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These  phenomena  can  be  especially  acute  in  the 
transonic-flow  regime,  where  shock-wave  positions 
and  strengths  are  sensitive  to  small  changes  in  the 
flow  parameters.  At  the  same  time,  the  development 
of  predictive  techniques  for  vortex  flows  in  this 
difficult  nonlinear  regime  has  not  kept  pace  with 
the  evolution  of  engineering  methods  for  flows  at 
either  higher  or  lower  Mach  numbers;  see,  for 
example,  the  recent  collection  of  papers  in  Ref.  1. 

A  particularly  complex  practical  example  of 
this  type  is  the  passage  of  helicopter  rotor  blades 
through  their  vortical  wakes  at  low  and  moderate 
flight  speeds.  The  blade  tips,  which  trail  strong, 
concentrated  tip  vortices,  trace  out  prolate 
cycloidal  paths  in  space,  leading  to  a  variety  of 
possible  blade-vortex  Interactions.  The  generic 
problem,  sketched  in  Fig.  1,  can  be  viewed  as  an 
unsteady,  three-dimensional  close  encounter  of  a 
curved-line  vortex,  at  an  arbitrary  Intersection 
angle  A,  with  a  high-aspect-ratio  lifting  surface 
that  is  executing  combined  rotational  and  transla¬ 
tional  motion  at  transonic  speeds.  The  limiting 
cases  for  A  -  0°  and  90°  are  illustrated  in 
Figs.  2a  and  2b,  respectively;  the  former  is 
essentially  two-dimensional  but  unsteady,  whereas 
the  latter  can  be  considered  as  steady  but  highly 
three-dimensional.  The  recent  paper  by  George  and 
Chang2  is  recommended  for  more  discussion  of  these 
representations  and  their  aeroacoustic  Implications. 

This  investigation  was  motivated  primarily  by 
the  situation  depicted  ir.  Fig.  2a,  as  a  first  step 
in  developing  practical  methods  for  predicting  gen¬ 
eral  three-dimensional,  unsteady  interactions.  In 
this  vein,  the  following  section  describes  the 
basic  formulation  of  the  problem  and  the  modifica¬ 
tion  of  an  existing,  efficient,  and  well-exercised, 
two-dimensional  computer  code3*4  to  treat  this 
class  of  unsteady,  transonic  vortex- interaction 
problems.  However,  the  manner  in  which  the  vortex 
is  introduced  readily  lends  itself  to  other  types 
of  convected  flow-field  disturbances,  such  as 
gusts,  and  this  extension  can  easily  be  combined 
with  the  inherent  capability  of  the  original  code 
to  treat  airfoils  with  moving  control  surfaces. 
Consequently,  the  method  has  also  been  applied  to 
cases  of  distributed  sinusoidal  gusts  in  the  free 
stream  and  to  the  alleviation  of  such  gusts  by 
means  of  control-surface  deflections,  in  addition 
to  the  concentrated  vortex-airfoil  Interaction 
problem.  Results  for  representative  examples  of 
each  of  these  problems  are  discussed  and  analyzed 
in  subsequent  sections  of  the  paper. 


II.  Numerical  Formulation 

The  transonic  small-disturbance  approximation 
to  the  velocity-potential  equation  was  chosen  for 
the  present  Investigation,  as  a  fruitful  compromise 
between  simplicity,  computational  efficiency,  and 
accuracy.  Viscous  effects  are  thus  ignored,  and 
the  results  are  also  acknowledged  to  be  inaccurate 
in  the  limstdiate  vicinity  of  the  leading  edge. 
However,  most  of  the  qualitative  features  of  the 
Interaction  phenomena  and  considerable  quantitative 
Information  can  be  obtained  within  the  scope  of 
this  approximation.  Also,  useful  solutions  can  be 
calculated  using  the  basic  algorithm  and  much  of 
the  original  coding  of  an  established  numerical 
program  called  LTRAN2  (Refs.  3  and  4).  However, 
the  special  features  of  the  vortex-and-gust- 
lnteractlon  problems  considered  in  this  paper  did 


require  several  additions  to  the  basic  equations 
and  boundary  conditions  and,  hence,  some  signifi¬ 
cant  modifications  to  the  code  LTRAN2,  as  outlined 
below. 


Governing 


auation  and  Numerical  Algorithm 


The  unsteady,  transonic  small-disturbance 
equation  for  the  velocity  potential  can  be  written 
in  the  following  form3*4: 


B»tt  +  2B*xt 


£  KCi  +  W»x)  + 


B  -  Mie"2'3 
C,  -  (1  -  M^e-2'3 
C2  -  -(Y  +  DM® 

V$  •  <J  -  1 

Here  e  is  the  airfoil  thickness  ratio,  and  the 
dimensionless  quantities  x,  y,  t,  <$,  and  $  are 
scaled  by  c,  ce_l  ,  c/U„,  U„,  and  ce2/sU„, 
respectively.  In  the  LTRAN2  codes,  the  exponent  m 
appearing  in  the  coefficient  C2  of  the  nonlinear 
term  $X4XX  can  be  set  to  the  value  2  for  the 
so-called  Spreiter,5  or  nonconservative,  scaling, 
or  allowed  to  vary  as  an  empirical  function  of 
local  and  free-stream  Mach  numbers.3  In  the  latter 
case,  Eq.  (1)  closely  approximates  the  formal 
small-disturbance  limit  of  the  conservative  form 
of  the  full-potential  equation. 6  Inclusion  of  the 
various  additional  nonlinear  terms,  for  example, 
*xx»  $x*xt>  which  arise  in  sundry  derivations  of 
the  small-disturbance  equations,  does  not  seem  to 
be  warranted  at  this  stage. 

On  the  other  hand,  the  term  B*tt,  which  was 
dropped  in  previous  versions  of  the  code  LTRAN2, 
should  be  included  here,  because  the  nondimenslonal 
time  scale  of  the  vortex  or  gust  Interaction, 

U„t/c,  may  be  of  the  order  of  unity.  This  has  been 
accomplished  by  adding  contributions  to  both 

the  x-  and  y-sweeps  of  the  alternating-direct ion- 
implicit  (ADI)  scheme  in  the  code  LTRAN2-HI 
(Ref.  4);  in  addition,  the  more  robust  monotone 
type-dependent  differencing  switch  of  Goorjlan  and 
Van  Buskirk7  was  substituted  for  the  standard 
Murman-Cole  switch  in  LTRAN2-HI.  Thus,  in  the 
present  case  the  solution  is  advanced  from  time- 
level  n  to  n  +  1  as  follows: 

x-sweep: 

(»*  _  2«n  +  *n_1)  + u*  _ 


V  +  6yy® 


y-sweep: 


— 2_  ($n+i  _  »*)  +  6  (♦n+1  -  **) 

(At)2  u  * 

■  i  V»n+1 '  *n)  <3) 

Here  $  is  understood  to  be  subscripted  with  i 
for  the  x-direction  and  j  for  the  y-direction. 


■ 


and  Is  Che  intermediate  value  of  the  unknown 

potential  at  the  end  of  the  predictor  step, 

Eq.  (2).  Following  approximately  the  notation  of 
Refs.  3  and  7, 

Dxf  -  tj  +  Axf 
?i-l/2,j  “  I  lC*“  +  (Ci  + 

^i-l/2, j  “  ~2  +  ^Ci  +  2c2u)tx**  j] 

t  r  »  lllihil  - 

^xr  *  1 

2  <*i+1  ’  Xi-!> 

t;  ^ i— i / 2 ,  j  ~  ^i— a/ 2 ,  j 

axr  *  1 

~2  ^xi+i  “  xl-i^ 


Vi.j 


~  *1-1, i  . 

xi  -  xi-i  ’ 
®in(u,  ui-i/j.j) 


Mi.j 


2  (xi+i  "  xl-i> 


u  *  max(u,  Ui_!/2tj) 

u  -  ~C1/C2 

Ui-X/2,j  ”  Vi.j 

5yy$i, j  *  2[<4i,j+i  ‘  'ti,j)/(yj+1  "  Fj>  -  ($i,j 

-  *itj-i)'Cyj  -  yj-xJi/Cyj+i  -  yj-^ 

It  should  be  mentioned  that  the  evaluation  of 


the  term 


in  this  new  scheme  is  only  first- 


order  accurate  in  both  time  and  space,  whereas  the 
basic  l. TRAN 2  algorithm  is  second-order  accurate  in 
time  and  first-order  accurate  in  space.  For  the 
applications  considered  in  this  paper,  however, 
this  reduction  in  theoretical  accuracy  is 
Inconsequential . 

Another  aspect  of  the  present  method  is  that 
by  Including  a  first  approximation  to  in  the 

x-sweep,  the  solution  at  time-level  n  -  1,  tn-1, 
does  not  appear  in  the  subsequent  y-sweep.  Con¬ 
sequently,  at  any  given  step  in  the  calculation, 
only  two  two-dimensional  arrays  are  required  in  the 
computer  memory  to  store  the  necessary  intermediate 
values  of  the  solution.  This  characteristic 
enabled  us  to  revise  the  LTRAN2-HI  code  to  solve 
Eqs.  (2)  and  (3)  without  having  to  store  any  addi¬ 
tional  large  data  arrays.  The  resultant  savings 
in  storage  requirements  can  be  an  important  advan¬ 
tage  for  computer  systems  with  limited  memory 
capacity,  and  it  would  be  even  more  so  in  three- 
dimensional  problems. 

Figure  3  shows  the  results  of  Including  the 
*tt  term  for  a  linear  test  case,  in  comparison 
with  the  results  obtained  with  the  original  code 
LTRAN2  (Ref.  3)  and  with  the  subsequent  refinements 
of  Hessenius  and  Goorjian"  in  LTRAN2-HI .  Particu¬ 
larly  noteworthy  is  the  improvement  in  the  pitching 
moment  results,  which  were  a  problem  before. 


B.  Introduction  of  Nonuniformities  in  the  Free 
Stream 

A  concentrated  potential  vortex  can  be  intro¬ 
duced  into  an  irrotational  flow  field  in  either  of 
two  different  ways,  as  illustrated  schematically 
in  Fig.  4.  The  first,  which  we  shall  call  the 
branch-cut  method,  is-  due  to  Caradonna6  and 
Caradonna  fit  al.®  and  is  a  straightforward  exten¬ 
sion  of  classical  potential  theory  for  lifting 
bodies.  As  indicated  in  Fig.  4a,  a  branch  cut  is 
introduce-'  between  the  vortex  and  the  outer  bound¬ 
ary  of  the  flow  field,  and  a  jump  in  potential 
equal  to  the  strength  of  the  vortex  is  prescribed 
across  this  branch  cut.  The  flow  remains  irrota¬ 
tional,  of  course,  outside  the  airfoil  and  vortex 
branch  cuts.  If  the  vortex  moves  through  the  flow 
field,  the  branch  cut  moves  accordingly. 

In  finite-difference  implementations  of  the 
branch-cut  method,  the  natural  choice  is  to  make 
the  vortex  and  airfoil  cuts  coincident  with  grid 
lines.  Then,  as  the  vortex  moves  through  the  flow 
field,  the  logic  of  the  code  must  move  the  branch 
cut  to  new  grid  lines;  but  this  must  be  done 
smoothly  to  avoid  introducing  artificial  distur¬ 
bances  into  the  solution  domain.  One  such  smooth¬ 
ing  method  was  recently  demonstrated  by  Murman  and 
Stremel,10  who  distributed  individual  point  vor¬ 
tices  over  neighboring  grid  points  using  Baker's 
vortex-in-cell  area-weighting  and  bilinear- 
interpolation  scheme.11 

The  Becond  approach,  which  we  adopt  here  and 
call  the  prescribed-disturbance  method,  evolved  in 
private  discussions  with  John  Steinhoff,  who 
recently  illustrated  the  method  in  Ref.  12  for  the 
steady,  three-dimensional,  full-potential  equation 
and  a  wing-vortex  interaction  problem.  The  basic 
idea,  indicated  schematically  in  Fig.  4b,  also 
parallels  certain  aspects  of  Goldstein’s13  and 
Kerschen  and  Myers'1"  aeroacoustic  extensions  of 
the  so-called  rapid-distortion  theory  of  turbu¬ 
lence,  using  the  linear  wave  equation.  Here  the 
essential  point  is  that,  although  Eq.  (1)  is  non¬ 
linear  and  Independent  solutions  are  not  super- 
posable,  the  velocity  field  can  still  be  split 
into  three  parts:  1)  the  (uniform)  free  stream, 

2)  a  prescribed  disturbance,  <fv,  for  example, 
gusts,  concentrated  vortices,  etc.,  and  3)  the 
disturbance  potential  Vif  due  to  the  airfoil 
itself.  That  is, 

+  v*  (4) 

If  Eq.  (4)  is  substituted  into  the  governing 
equation(s),  a  new  and  more  complicated  equation 
will  result,  unless  the  original  equation  was 
linear.  Moreover,  the  prescribed  disturbance  at 
each  time-step  or  iteration  may  turn  out  to  depend 
on  the  solution  of  the  governing  equation  at  the 
previous  time  step  or  iteration.  However,  it  may 
also  turn  out  that  a  given  solution  technique  or 
algorithm  that  has  been  developed  for  the  original 
governing  equation  will  work  equally  well  for  the 
new,  more  complicated  equation.  That  is  the  case 
here,  which  is  particularly  straightforward  in  the 
small-disturbance  limit. 

For  convenience  and  for  purposes  of  illustra¬ 
tion,  “  1^  +  }vv  is  chosen  to  satisfy  the 
linearized  version  of  Eq.  (1)  (i.e.,  with  C2  ■  0) . 
Then  substituting  Eq.  (4)  into  Eq.  (1)  and 


3 


subtracting  out  the  linear  terms  gives  the  follow¬ 
ing  equation  for  the  unknown  disturbance  potential 

B^t  +  2B*xt  -  +  C,  ^  +  Uy)2  +  ♦yy 

(5) 

The  same  substitution  is  also  made  in  the 
original  boundary  conditions,  which  were  1)  flow 
tangency  on  the  body,  2)  no  disturbances  at 
x  ■*  and  y  -*■  ±®,  and  3)  Cp  ■  0  at  x  +■>. 

The  new  small-disturbance  boundary  condition  on  the 
body,  yj,  ■  F(x,t),  becomes 

*y(x,y  -*■  0)  -  Fx  +  Ft  -  vv  (6) 

The  wake  condition  of  zero  pressure  jump  across  the 
airfoil  branch-cut  remains  unchanged,  that  is, 

(rx  +  rt)w  -  o  (?) 

With  the  respecification  of  the  appropriate 
boundary  conditions  at  the  outer  boundaries  of  the 
computational  domain,  Eqs.  (5) — ( 7)  result  in  a  new 
well-posed  problem  for  $  that  can  be  solved  by 
the  same  basic  algorithm  as  before.  Equations  (2) 
and  (3)  read  the  same,  except  it  is  understood  that 
$  is  replaced  by  $,  the  various  u,  u,  u,  etc. 
are  replaced  everywhere  by  (u  +  uv) ,  etc.,  and  the 
nonlinear  term  becomes 

Dxf  "  ^x*i-l/2,J  +  ^x*i-i/2,j  +  "2  C2l5x[(5x*n  +  “^“vl 

(8) 


small-disturbance  Bernoulli  equation  for  Cp  Is 
second  order,  and  hence  negligible. 

In  contrast  with  the  branch-cut  method,  we 
see  that  the  present  approach  introduces  additional 
nonlinear  terms  into  the  governing  equation  and 
into  the  boundary  conditions  on  the  airfoil,  but 
that  the  moving  branch-cut  logic  of  the  former 
method  is  effectively  replaced  by  the  simple 
expression  of  Eq.  (lib).  In  addition,  the  pre¬ 
scribed  disturbance  need  not  be  restricted  to  the 
form  of  Eqs.  (9)  and  (10);  for  example,  a  vortex 
with  a  viscous  core  and  finite  velocity  at  its 
origin  could  also  be  specified.  In  fact,  this  was 
found  to  be  necessary  for  most  of  the  examples 
considered,  as  explained  in  Sec.  IIC. 

The  prescribed-disturbance  method  can  also  be 
extended  to  much  more  general  cases  of  vortical  or 
other  disturbances.  This  is  especially  simple  in 
the  small-disturbance  limit,  where  most,  or  per¬ 
haps  all,  of  the  resulting  cross  terms  will  be 
negligible  and  the  coding  modifications  will  be 
minimal.  Typical  practical  examples  include  small- 
amplitude  atmospheric  gusts,  of  either  the  sharp- 
edged  or  traveling-wave  type,  which  heretofore  have 
only  been  treated  with  linear  analyses.  Further¬ 
more,  if  the  gust  is  uniform  in  the  direction 
normal  to  its  disturbance  velocity,  its  effect 
enters  only  through  the  airfoil  boundary  condition, 
Eq.  (6),  if  it  is  a  vertical  gust,  or  through  Cp 
and  the  nonlinear  term  (d>x  +  uq)2  in  Eq.  (5)  if  it 
is  a  horizontal  gust.  Examples  of  the  former  type 
are  given  in  Sec.  IV,  with  and  without  accompanying 
control-surface  deflections. 

C.  Effects  of  Secondary  Numerical  Parameters 


In  this  derivation,  the  strength  and  structure 
of  the  prescribed  disturbance  are  implicitly 
assumed  to  be  unaltered  by  the  presence  of  the 
airfoil,  and  it  is  convected  through  the  flow 
field  at  the  velocity  of  the  local  fluid  elements. 
Then,  for  the  particular  case  in  which  the  imposed 
disturbance  is  a  concentrated,  moving  potential 
vortex,  its  disturbance  field  and  position  in 
scaled  variables  are  given  by 


Calculations  were  performed  with  both 
m  «  f(M,M„)  and  with  m  *  2  in  Eq.  (1),  and  with 
certain  other  parameters  in  the  numerical  code 
that  were  chosen  somewhat  arbitrarily.  These 
Included  the  step  size  At;  the  number  of  grid 
points  and  their  distribution;  and  the  value  of 
x0  where  the  vortex  was  initially  introduced  and 
the  radius  a  of  the  core  of  the  vortex,  in  the 
airfoil-vortex  interaction  examples. 


“v 


JLf _ 

2™  |_(x  -  Xl 


r  +  (y 


yi>2e 


-2/s] 


(9) 


r  r _ y  -  Ji _ 

V  2lre  |_(x  -  Xj)*  +  (y  -  yj)2e_2/3 


(10) 


The  effect  of  the  exponent  m  in  Eq.  (1)  was 
found  to  be  merely  what  would  be  expected  from  the 
steady  case  for  Spreiter  scaling  versus  fully  con¬ 
servative  scaling.  That  is,  the  former  scaling 
tended  to  produce  weaker  shock  waves  and  smaller 
supersonic  zones.  However,  the  basic  qualitative 
effects  of  the  airfoil-vortex  or  airfoil-gust 
interactions  were  unaffected. 


n 

3  +  £  <Si_  At  (lib) 

m“i 

where  •  (xI,y1)  is  the  instantaneous  position 
of  the  vortex;  is  the  original  position  of  the 
vortex  when  the  calculation  is  initiated;  T  is 
the  vortex  strength,  normalized  by  c  and  Ua;  and 
*  (uj.Vj)  is  the  velocity  at  the  vortex  location 
tlt  obtained  from  the  numerical  solution  of  Eq.  (5) 
at  time-level  m.  We  note  in  passing  that  although 
such  a  vortex  contributes  first-order  terms  in 
Eqs.  (5)  and  (6) ,  its  net  contribution  to  the 


Most  of  the  calculations  performed  in  this 
investigation  were  done  with  the  standard  LTRAN2 
grid  of  113  points  in  the  x-direction  and  97  in 
the  y-direction,  with  smooth  stretching  in  both 
directions  to  outer  boundaries  200  chord  lengths 
away.  However,  finer  x-grids  of  up  to  186  points 
were  found  to  be  necessary  for  adequate  resolution 
of  the  pressure-wave  propagation  in  some  of  the 
airfoil-vortex  interaction  cases  with  the 
NACA  64A006  airfoil.  These  finer  grids  were  also 
required  for  the  oscillating-airfoil  and  sinusoidal 
gust  cases  at  reduced  frequencies  greater  than 
about  two,  where  the  computed  results  began  to 
deviate  from  linear  theory  otherwise.  The  compu¬ 
tation  times  were  dependent  on  the  number  of  grid 
points;  CPU  times  were  approximately  1.5  *  10“5  sec 
per  grid  point  per  time  step  on  the  Cray  1-S 
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computer  at  Ames  Research  Center  (or  between  1  and 
5  min  per  case) ,  and  about  twice  this  on  the 
CDC  7600  machine. 


The  calculations  of  pitching  oscillations  and 
gust  responses  were  normally  done  with  240  time 
steps  per  cycle  over  5  or  6  cycles  of  oscillation, 
after  it  was  established  that  the  solutions  were 
independent  of  At  for  more  than  abouc  180  steps 
per  cycle.  Most  of  the  airfoil-vortex  interaction 
cases  were  performed  with  At  «  0.02.  The  accuracy 
of  the  nonlinear  results  began  to  degrade  for 
At  «  0.05,  and  spurious  oscillations  in  the  solu¬ 
tions  sometimes  appeared  for  At  <  0.02. 


The  results  in  Sec.  Ill  show  that  the  vortex 
began  to  influence  the  pressure  distribution  on  the 
airfoil,  although  weakly,  when  it  was  very  far 
upstream.  For  computational  economy,  therefore, 
the  vortex  was  normally  introduced  at  7  to  10  chord 
lengths  upstream  of  the  leading  edge,  and  then  its 
strength  was  increased  linearly  with  time  from 
zero  initially  to  the  full  value  after  1  chord 
length  of  travel.  Solutions  obtained  with 
x0  -  -7,  -10,  and  -20  were  found  to  be  virtually 
the  same  by  the  time  the  vortex  was  within  3  chord 
lengths  of  the  leading  edge  and  the  pressures  on 
the  airfoil  began  to  vary  rapidly. 


Numerical  instability  problems  were  encoun¬ 
tered  early  in  the  investigation  when  the  potential 
vortex  of  Eqs.  (9)  and  (10)  was  employed  in  non¬ 
linear  calculations.  These  difficulties  were 
traced  to  excessive  values  of  Uy  and  its 
.-derivative;  consequently,  a  vortex  with  a  finite 
core  velocity  was  introduced,  having  an  induced 
velocity  field  given  by 


r 

2irer 


(12) 


where 


introducing  the  vortex,  the  small-disturbance  equa¬ 
tions  with  Spreiter  scaling,  and  a  somewhat  differ¬ 
ent  ADI  algorithm.  Accordingly,  the  present  code 
was  first  checked  out  for  the  cases  published  in 
Ref.  9,  for  the  NACA  0012  airfoil  at  M„  =  0.8. 

Some  representative  examples  are  given  in  Figs.  6 
and  7.  In  this  particular  case,  the  vortex  is 
prescribed  to  move  at  free-stream  velocity  along 
the  path  yj  *  constant,  and  Xj  =  xt  +  UTOt/c, 
following  Ref.  9.  The  effect  of  this  approxima¬ 
tion  is  examined  in  Sec.  IIIC. 

Figure  6  shows  the  pressure  distribution  on 
the  airfoil  at  several  stages  of  the  vortex  pas¬ 
sage.  The  agreement  between  the  two  methods  is 
satisfactory,  except  for  the  time  when  the  vortex 
is  beneath  the  midchord  of  the  airfoil.  It  should 
be  mentioned  that  for  this  particular  set  of  con¬ 
ditions,  the  interaction  is  quite  severe  in  terms 
of  the  small-disturbance  approximation.  The 
numerical  solutions  seem  to  be  marginally  stable, 
depending  on  the  mesh  distribution,  on  At,  and  on 
the  vortex  core  radius.  Nevertheless,  the  two 
independent  sets  of  results  illustrate  most  of  the 
essential  features  of  the  airfoil-vortex 
interaction. 

In  a  very  general  sense,  the  initial  develop¬ 
ment  of  the  pressure  difference  across  the  airfoil 
resembles  the  expected  response  of  the  flow  to  a 
decrease  in  angle  of  attack;  that  is,  the  vortex 
induces  a  spacially  varying  "downwash,"  or  negative 
v-velocity  perturbation,  in  the  flow  approaching 
the  airfoil.  However,  the  effect  is  much  larger 
on  the  lower  surface  than  on  the  upper,  a  result 
of  nonlinear  effects.  Also,  the  actual  shape  of 
the  pressure  distributions,  the  shock-wave 
strengths  and  positions,  and  the  time  history  of 
the  airfoil  pressures  are  quite  different,  even 
qualitatively,  from  the  results  for  airfoil  oscil¬ 
lations,  or  from  the  results  for  the  sinusoidal 
gusts  considered  in  Sec.  IV. 


r2  -  (x  -  xx)2  +  (y  -  yi)2e-2/3 

This  velocity  field  is  illustrated  in  Fig.  5.  It 
is  noteworthy  that  qv  is  very  nearly  equal  to  the 
ideal-vortex  value  of  Eqs.  (9)  and  (10)  for  r  >  2a. 
Therefore,  the  predominant  effect  of  a  is  felt 
only  very  locally  in  the  flow  field,  and  there  only 
through  the  nonlinear  term  ($x  +  Uy) 2 ,  in  the  small- 
disturbance  limit  of  the  present  study.  Test  cal¬ 
culations  showed  the  solutions  to  be  essentially 
independent  of  a  for  the  range  0.02  <  a  <  0.07; 
the  value  0.05  was  used  for  the  results  presented 
in  Sec.  III.  The  resolution  of  this  numerical 
difficulty  by  means  of  a  vortex  with  a  core  may 
have  implications  for  future  work  using  the  branch- 
cut  method,  as  well. 


III.  Airfoil-Vortex  Interactions 


Finally,  after  the  vortex  passes  behind  the 
airfoil,  its  induced  field  becomes  an  "upwash" 
distribution  of  vertical  velocity.  The  qualita¬ 
tive  effect  of  this  is,  very  approximately,  the 
inverse  of  the  behavior  when  the  vortex  is  upstream 
of  the  airfoil.  However,  the  return  of  the  flow 
around  the  airfoil  to  the  original  state  is 
extremely  slow,  even  though  the  lingering  effects 
are  much  weaker. 

B.  Results  for  the  NACA  64A006  Airfoil 

This  airfoil,  which  has  been  used  in  numerous 
numerical  and  experimental  studies,  is  not  only 
thinner  than  the  NACA  0012  section,  but  has  a  sig¬ 
nificantly  smaller  leading-edge  radius  as  well. 
Therefore,  it  might  be  expected  to  be  more  sensi¬ 
tive  to  the  vortex-induced  downwash  than  the 
NACA  0012  profile. 


A.  Results  for  the  NACA  0012  Airfoil 

The  problem  of  a  concentrated  potential  vortex 
convectlng  past  a  stationary  airfoil  is  indicated 
schematically  in  Fig.  2a.  Parathasarathy1 5  and 
Chow  and  Huang2 *  have  found  incompressible  solu¬ 
tions  using  conformal  mapping  techniques  for  this 
general  case.  As  discussed  in  Sec.  II,  this  prob¬ 
lem  was  first  treated  in  the  transonic  case  by 
Caradonna  et  al.,9  using  the  branch-cut  method  for 


Figures  8-11  show  the  results  for  an  unsteady 
case  with  m  «  2  and  with  the  same  vortex  strength 
and  initial  value  of  yv^  as  in  the  previous 

example,  but  at  M„  *  0.85  and  with  a  force-free 
trajectory  given  by  Eq.  (11).  Here  the  airfoil 
Incidence  is  zero,  and  the  dramatic  differences 
between  Cp  on  the  upper  and  lower  surface  are 
due  solely  to  the  vortex  interaction.  As  before, 
the  vortex  distorts  the  flow  on  the  lower  surface 
of  the  airfoil  more  than  on  the  upper. 


5 


As  In  the  previous  case,  the  effect  of  the 
vortex  on  the  airfoil  flow  field  is  felt  for  vortex 
positions  many  chord  lengths  upstream  and  down¬ 
stream  of  the  origin.  However,  the  most  rapid  and 
dramatic  changes  in  the  airfoil  pressure  distribu¬ 
tions  occur  when  the  vortex  is  within  1  or  2  chord 
lengths  of  the  airfoil.  For  this  particular  air¬ 
foil,  these  principal  variations  in  Cp  have  two 
other  distinct  stages  of  pressure  fluctuations 
superimposed  upon  them.  The  first  is  the  rapid 
development  of  a  large  suction  peak  near  the  lower 
leading  edge  as  the  vortex  approaches.  This  is 
followed  soon  thereafter  by  an  even  more  rapid 
increase  in  Cp,  that  is,  by  a  rapid  compression. 

As  noted  by  George  and  Chang,2  a  compression  wave 
seems  to  leave  the  airfoil  at  this  point  and  propa¬ 
gate  into  the  oncoming  flow.  The  inverse  of  this 
phenomenon  occurs  in  the  upper-surface  leading-edge 
region,  that  is,  a  local  compression  followed  by  an 
expansion.  However,  the  magnitude  and  rate  of 
change  of  the  fluctuating  pressure  there  is  con¬ 
siderably  less,  and  an  expansion  wave  does  not  seem 
to  propagate  forward. 

The  second  pattern  of  moving  pressure  waves 
occurs  some  time  after  the  vortex  passes  downstream 
of  the  trailing  edge.  The  shock  wave  on  the  lower 
surface  begins  to  creep  forward,  weakening  as  it 
moves.  Eventually  a  weak  compression  wave  moves 
forward,  ahead  of  the  shock  wave,  but  it  essen¬ 
tially  vanishes  before  it  reaches  the  leading  edge. 
This  lower  surface  behavior  resembles  the  so-called 
type  B  unsteady  shock-wave  motion  of  Tijdeman,17  as 
noted  by  George  and  Chang.2  During  this  time,  the 
upper-surface  shock  wave  gradually  weakens,  as  the 
flow  returns  slowly  toward  the  initial  undisturbed 
state.  Although  they  move  significantly  and  change 
in  strength,  the  primary  shock  waves  on  both  the 
upper  and  lower  surfaces  remain  attached  to  their 
respective  sides  of  the  airfoil  throughout  the 
vortex-encounter  events. 

Figures  12  and  13  show  the  same  type  of 
results  as  Figs.  8  and  9,  but  for  a  lower  free- 
stream  Mach  number.  This  is  the  airfoil  and  Mach 
number  for  which  the  Tijdeman  type  C  shock-wave 
motion  was  calculated  by  Ballhaus  and  Goorjian3 
for  an  oscillating  flap.  A  weak  type  C  behavior 
on  the  lower  surface  is  evident  here  as  well,  but 
the  compression  wave  essentially  vanishes  before 
reaching  the  leading  edge. 

C.  Effects  of  Various  Approximations 

A  number  of  test  cases  were  run  for  the  pre¬ 
ceding  flow  conditions,  with  various  terms  deleted 
one  at  a  time  from  the  basic  equations.  The 
results  of  this  numerical  experimentation  are  sum¬ 
marized  in  Figs.  14-19,  which  show  the  time-history 
of  the  lift  and  two  instantaneous  pressure  distri¬ 
butions  for  each  caBe. 

Figure  14  shows  the  effect  of  neglecting  the 
4tt  term;  as  reported  by  George  and  Chang,2  this 
term  does  not  seem  to  be  as  important  as  was 
originally  thought.  However,  it  was  found  to  have 
somewhat  of  a  stabilizing  effect  on  calculations 
that  were  marginally  stable. 

Figures  15  and  16  show  the  effect  of  neglect¬ 
ing  the  horizontal  and  vertical  components,  respec¬ 
tively,  of  the  vortex-induced  velocity.  It  is 
clear  from  these  two  figures  that  the  latter  is 
the  most  important,  at  least  for  this  particular 


Mach  number  and  airfoil  section.  This  result 
indicates  that  the  approximation  of  George  and 
Chang2  in  neglecting  Uy  was  not  too  serious. 
However,  in  comparison  with  the  results  in 
Figs.  8-11,  the  perturbations  in  the  lower-surface 
pressures  in  the  leading-edge  region  are  too  large 
when  uv  is  neglected. 

Figure  17  shows  the  effect  of  specifying  that 
the  vortex  move  along  a  straight  line,  yj  «  y0 
and  xx  ■  x0  +  Umt/c,  rather  than  allowing  it  to 
move  on  the  force-free  path  given  by  Eq.  (11).  In 
the  latter  instance,  the  vortex  is  displaced  far¬ 
ther  from  the  airfoil,  and  its  effect  on  the  flow 
field  around  the  airfoil  weakens  as  a  result. 

Finally,  Figs.  18  and  19  show  the  effect  of 
neglecting  the  nonlinear  term  altogether,  that  is, 
a  linear  calculation.  The  overall  lift  and  moment 
are  not  too  different  from  the  values  calculated 
using  the  full  nonlinear  formulation,  but  the  time 
history  of  the  leading-edge  pressures  shows  much 
less  rapid  fluctuations.  The  changes  in  the  upper 
and  lower  pressures  from  the  initial  values  with¬ 
out  the  vortex  present  are,  of  course,  exactly 
equal  and  opposite  for  this  case. 


IV.  Airfoil-Gust  Interactions 

The  modifications  outlined  in  Sec.  II  make  it 
possible  to  calculate  the  effects  of  gusts  in  the 
free  stream.  This  is  illustrated  for  the  particu¬ 
lar  case  of  the  classical  longitudinally  traveling, 
sinusoidal  gust  of  vertical  velocity, 

VG  “  wi  sinjut  -  kg(x  -  1/2)]  (13) 

This  expression  for  vG  replaces  vv  in  the 
boundary  condition  on  the  airfoil,  Eq.  (6).  There 
is  no  u-component  of  velocity  for  this  particular 
type  of  gust.  Also,  the  distortion  of  the  gust 
field  by  the  presence  of  the  airfoil,  as  treated 
by  Kerschen  and  Myers11*  in  the  linear  case,  is 
ignored . 

Figure  20  shows  linear  results  obtained  with 
the  present  code  for  the  Sears  lift  function18  in 
comparison  with  the  compressible  results  of 
Graham19  for  low-  and  high-subsonic  Mach  numbers, 
over  a  wide  range  of  reduced  frequencies.  As  in 
Fig.  3,  k  ■  wc/U«,  is  based  on  the  chord  length, 
and  here  kg  is  inversely  proportional  to  the 
wavelength  of  the  gust.  The  agreement  is  com¬ 
parable  to  that  shown  in  Fig.  3  for  pitching 
oscillations. 

A.  Nonlinear  Gust  Response 

Figure  21  shows  the  distortions  in  the  pres¬ 
sure  distribution  on  an  NACA  0012  airfoil  in  tran¬ 
sonic  flow  due  to  a  sinusoidal  gust.  In  this 
case,  the  fluctuations  in  the  solution  are  anti¬ 
symmetric  on  the  upper  and  lower  surfaces  of  the 
airfoil.  The  steady  and  first-harmonic  content  of 
these  fluctuations  are  compared  with  linear 
results  in  Figs.  22  and  23  at  M«,  ■  0.80  and 
k  •  0.2  and  2.0,  respectively.  The  peak  in  the 
nonlinear  results  near  midchord  is  a  direct  result 
of  the  motion  of  the  shock  wave;  see,  for  Instance, 
Ref.  20.  This  effect  attenuates  with  increasing 
reduced  frequency,  similar  to  the  trends  of  pitch¬ 
ing  or  plunging  airfoil  oscillations. 
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In  various  cases  summarized  In  Figs.  20-23. 

1c  Is  seen  that  the  magnitudes  of  the  lift  and 
pitching  moment  decrease  significantly  and  mono- 
tonically  with  increasing  gust  frequency.  This 
might  be  expected  from  the  well-known  trends  of  the 
linear  case,  but  it  is  very  different  from  the 
effect  of  oscillatory  airfoil  motion.  For  pitching 
or  plunging  oscillations  in  transonic  flow,  the  air 
loads  typically  decrease  in  magnitude  with  increas¬ 
ing  frequency  at  low  values  of  the  reduced  fre¬ 
quency  k,  but  increase  with  frequency  at  higher 
values  of  k.  However,  the  limiting  low-frequency 
behavior  of  all  three  becomes  the  same  as  k  -*■  0. 

B.  Combined  Gusts  and  Control-Surface  Motions 

A  number  of  calculations  were  performed  to 
examine  the  possibility  of  canceling  gust-induced 
unsteady  air  loads  by  the  appropriate  motion  of  a 
simple  trailing-edge  flap,  using  the  combined  capa¬ 
bility  of  the  original  LTRAN2  code  and  the  present 
modifications  for  calculating  gust  response.  The 
procedure  is  indicated  schematically  in  Fig.  24. 

For  example,  if  the  problem  were  linear  and  the 
gust  were  described  by  Eq.  (13),  then  the  fluctuat¬ 
ing  lift  due  to  the  gust  could  be  canceled  by 
deflecting  the  flap  as  follows: 

«F  -  A  sin(t  -  6)  (14) 

where  A  ■  G/F;  0  »  0G  -  0F;  G  and  0G  are  the  mag¬ 
nitude  and  phase  of  the  unsteady  lift  due  to  the 
gust  but  without  flap  deflections;  and  F  and  6p 
are  the  magnitude  and  phase  of  the  unsteady  lift 
caused  by  a  small  flap  deflection  (e.g.,  1°),  and 
reduced  frequency  kG  in  a  uniform  oncoming  flow. 

In  other  words,  the  principle  of  linear  superposi¬ 
tion  could  be  used  to  determine  the  requisite  flap 
deflection. 

In  the  nonlinear,  transonic  case,  however,  the 
foregoing  strategy  would  not  be  expected  to  work 
unless  the  gust  were  weak  and  the  shock-wave  motion 
in  response  to  the  gust  were  negligible.  The  cases 
depicted  in  Figs.  20-22,  for  example,  do  not  satisfy 
these  requirements.  Consequently,  when  Eq.  (14)  was 
implemented  for  these  conditions  of  and  kc,  sub¬ 
stantial  unsteady  air  loads  remained,  both  Cl  and  Cm- 

However,  it  was  found  that  one  additional 
iteration  of  the  preceding  strategy  was  adequate, 
that  is,  by  employing  the  following  flap  deflection: 

6p  «  Aj  sin(t  +  0j)  +  A2  sin(t  +  02)  (15) 


®G,  ~  ®F 


First,  a  calculation  was  done  using  Eq. 


specify  the  amplitude 


and  phase 


(14)  to 
.  Then  the 


resulting  response,  G2  and  0G^  were  used  in  a 

second  calculation  to  determine  A2  •  G2/F  and 
®2  "  ®G  “  ®F*  For  an  example  corresponding  to  the 

conditions  of  Fig.  21,  this  procedure  reduced  the 
magnitude  of  the  unsteady  lift  coefficient  from 
0.080  with  no  flap  deflection  to  0.002.  The 
resultant  pressure  distributions  are  shown  in 


25  in  harmonic  format.  In  particular,  the 


in-phase  component  of 


shows  one  spike  centered 


on  the  mean  shock-wave  position  and  another  near 
the  hinge  point  of  the  flap. 

Figure  25  also  shows  that  even  though  the  lift 
was  essentially  canceled,  the  pressure  fluctuations 
remained  large;  and,  in  fact,  a  significantly 
larger  pitching  moment  was  developed  with  flap 
oscillation  than  without.  Conversely,  applying 
the  strategy  of  Eq.  (15)  to  minimize  the  pitching 
moments  virtually  canceled  C^,  but  at  the  expense 
of  large  fluctuating  lift.  A  similar  perversity, 
but  of  a  somewhat  lesser  degree,  was  found  by 
Ballhaus  and  Goorjian21  in  their  attempts  to  cancel 
the  effects  of  plunging  motion  with  flap  deflec¬ 
tion;  they  showed  that  a  combination  of  leading- 
and  trailing-edge  flaps  was  required  to  cancel  both 
Cl  and  Cjj.  The  basic  problem  in  both  cases  is  that 
the  trailing-edge  flap  deflection  produces  an  aft 
loading  on  the  airfoil,  whereas  the  loading  result¬ 
ing  from  a  gust  or  from  plunging  oscillations  tends 
to  be  effectively  centered  much  closer  to  the 
quarter  chord,  x  =  0.25.  Therefore,  it  seems  to 
be  inherently  impossible  to  provide  effective 
active  control  of  gust  response  with  a  simple 
trailing-edge  flap.  On  the  other  hand,  the  pres¬ 
ent  computational  capability  provides  a  new  tool 
with  which  other  gust-alleviation  schemes  can  be 
studied. 


and  Conclusions 


Relatively  straightforward  modifications  and 
improvements  have  been  made  in  the  computer  code 
LTRAN2  that  permit  new  classes  of  unsteady 
transonic-flow  problems  to  be  solved  efficiently. 
Test  cases  have  validated  the  new  code,  and 
sample  results  illustrate  the  essential  effects 
of  convected  disturbances,  both  distributed  and 
concentrated,  on  the  flow  in  the  vicinity  of  an 
airfoil. 

The  calculated  results  for  discrete  vortex- 
airfoil  interactions  show  major  distortions  in  the 
flow  field  as  the  vortex  passes  by.  The  most 
important  factor  in  producing  these  distortions 
seems  to  be  the  time-  and  space-dependent  vertical 
velocity,  or  local  downwash  and  upwash,  that  is 
induced  by  the  moving  vortex,  although  its  hori¬ 
zontal  induced-velocity  component  also  contributes 
in  transonic  cases.  For  vortex  trajectories  pass¬ 
ing  close  to  the  airfoil,  it  is  important  to  allow 
the  vortex  to  move  along  a  force-free  path. 

The  time  histories  of  the  airloads  that  were 
obtained  by  running  the  code  in  the  linear  mode 
were  somewhat  similar,  qualitatively,  to  the  tran¬ 
sonic  results,  even  though  the  instantaneous  chord- 
wise  pressure  distributions  were  significantly 
different.  Also,  the  maximum  rates  of  change  of 
the  linear  pressure  results  were  significantly 
less. 

In  the  case  of  a  thin  profile  with  a  small 
leading-edge  radius,  large  and  rapidly  varying 
pressure  fluctuations  were  found  to  occur  near  the 
lower-surface  leading  edge.  Since  this  is  the 
region  in  which  the  small-disturbance  approxima¬ 
tion  breaks  down,  the  present  results  should  be 
checked  with  a  more  accurate  formulation. 
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The  capability  of  the  code  to  treat  distributed 
vortical  disturbances  in  the  free  stream  was  demon¬ 
strated  for  the  case  of  a  sinusoidal  gust  of  ver¬ 
tical  velocity.  Again,  the  trends  of  the  transonic 
results  were  qualitatively  similar  to  the  predic¬ 
tions  of  linear  theory,  but  important  quantitative 
differences  were  found,  as  would  be  expected,  when 
the  shock  waves  on  the  airfoils  moved  fore  and  aft. 

Combined  gust  response  and  oscillating  control 
surfaces  were  also  examined,  and  a  simple  strategy 
was  developed  for  suppressing  the  gust-induced  air 
loads  by  activating  a  trailing-edge  flap.  However, 
the  chordwise  pressure  distributions  associated 
with  gusts  are  very  different  from  those  produced 
by  trailing-edge  flap  oscillations.  Consequently, 
this  simple  approach  cannot  cancel  both  the  gust- 
induced  fluctuating  lift  and  the  pitching  moments, 
even  though  one  or  the  other  can  be  suppressed 
almost  completely  by  oscillating  the  flap  with  the 
proper  phase  and  amplitude. 
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Fig.  6 
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Instantaneous  pressure  distributions  during  an  airfoil-vortex  interaction:  NACA  0012  airfoil 
M«,  -  0.80,  a.  »  0.5°,  V  -  0.20,  yv  -  -0.26. 


VORTEX  POSITION.  xy 

Fig.  7  Lift  and  pitching-moment  variations  versus 
instantaneous  x-position  of  the  vortex;  same  con¬ 
ditions  as  in  Fig.  6. 
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Fig.  12  Instantaneous  pressure  distributions  during  an  airfoil-vortex  interaction:  NACA  64A006  airfoil 

H»  “  0.822,  a  -  0,  T  =  0.20,  yVo  -  -0.26. 


Fig.  13  Lift,  pitching  maaent,  and  vortex  dis- 
placement  for  the  NACA  64A006  airfoil  at 

m_  -  n.H??. 
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19  Linear  calculation  of  surface  pressures. 


Fig.  20  Linear  calculation  of  the  lift  response 
of  an  airfoil  to  a  sinusoidal  gust. 
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Fig.  21  Instantaneous  pressure  distributions  due  to  a  sinusoidal  gust:  NACA  0012  airfoil,  M„ 

a  =  0,  kG  -  0.50,  wj  =  1°. 
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